SCATTERING THEORY OF DISCRETE (PSEUDO) LAPLACIANS 

ON A WEYL CHAMBER 
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Abstract. To a crystallographic root system we associate a system of multi- 
variate orthogonal polynomials diagonalizing an integrable system of discrete 
pseudo Laplacians on the Weyl chamber. We develop the time-dependent 
scattering theory for these discrete pseudo Laplacians and determine the corre- 
sponding wave operators and scattering operators in closed form. As an appli- 
cation, we describe the scattering behavior of certain hyperbolic Ruijsenaars- 
Schneider type lattice Calogero-Moser models associated with the Macdonald 
polynomials. 
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1. Introduction 

A fundamental property of the solitonic solutions of integrable nonlinear wave 
equations is that their multi-particle scattering process decomposes into pairwise 
two-particle interactions |SCMI IXHl IN-ZI INI IFT| . This phenomenon is preserved at 
the quantum level: the corresponding solitonic quantum field theories are charac- 
terized by an iV-particle scattering matrix that factorizes in terms of two-particle 
scattering matrices [MJ IKBI) . As it turns out, this type of factorization can be 
understood heuristically as being a consequence of the integrability of the models 
in question |Ku[IRSc) . 

An archetype example of an integrable system with factorized scattering is the 
celebrated nonlinear Schrodinger equation (NLS). The quantum version of this 
model boils down to a bosonic A'^-particle system with a pairwise interaction via 
delta-functional potentials. The factorization of the scattering manifests itself 
through the asymptotics of the wave function, which is characterized by (prod- 
ucts of) two-particle scattering matrices (or c-functions) [CTlKlallTTxllKBII . 

In recent work, Ruijsenaars constructed a remarkably large class of quantum 
integrable lattice models of iV-particles exhibiting factorized scattering IR4j . The 
discrete systems in question arise by interpreting recurrence relations (or Fieri for- 
mulas) for symmetric multivariate orthogonal polynomials as quantum eigenvalue 
equations. Here the polynomial variable plays the role of the spectral parame- 
ter and the index (i.e. partition) labelling the polynomials is thought of as the 
discrete spatial variable. By analyzing the asymptotics of the polynomials as the 
degree tends to infinity, Ruijsenaars demonstrated that — for factorized orthogo- 
nality measures subject to certain technical conditions ensuring that the particle 
interaction is short-range and the spectrum is absolutely continuous — the corre- 
sponding discrete models are governed by a scattering matrix that factorizes into 
two-particle scattering matrices. An interesting particular case is that of the Mac- 
donald polynomials |M2| . The corresponding iV-particle model can be identified as 
a hyperbolic Ruijsenaars-Schneider type lattice Calogero-Moser system |R1I IR3j . 
At the level of classical Hamiltonian mechanics, the scattering of the corresponding 
integrable system was studied in great detail in Ref. )B,2| . 

It is known that the root systems of simple Lie algebras form a fruitful context for 
understanding Calogero-Moser systems and particle models with delta-functional 
potentials |OP[ IGul IHSL IHOL lU) . From this perspective, it is natural to ask for a 
generalization of Ruijsenaars' construction to the case of arbitrary root systems. 
The purpose of the present paper is to provide such a construction. 

More specifically, we associate to a crystallographic root system a system of 
Weyl-group invariant multivariate orthogonal polynomials on the Weyl alcove, char- 
acterized by a weight function that factorizes over the roots (of the root system) 
in terms of one-dimensional c-functions. The orthogonality implies that the poly- 
nomials satisfy a system of recurrence relations (Fieri formulas) . These recurrence 
relations are interpreted as eigenvalue equations for an integrable system of discrete 
pseudo Laplacians on the Weyl chamber. We develop the time-dependent scatter- 
ing theory for these discrete pseudo Laplacians and determine the corresponding 
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wave operators and scattering operators in closed form. For a specific choice of the 
weight function, our polynomials amount to the Macdonald polynomials associated 
with root systems iM3l IM4j . Again the corresponding integrable lattice model then 
permits identification as a discrete hyperbolic Ruijsenaars-Schneider type Calogero- 
Moser system jRll IR3I IDlj . For the type A root systems the Weyl group is the 
symmetric group and we reproduce the results of Ruijsenaars |H 4| . 

The wave- and scattering operators computed in this paper compare the dynam- 
ics generated by the discrete pseudo Laplacian to that of a free discrete Laplacian 
(corresponding to the case that the c- functions reduce to constant functions). Our 
study of the scattering consists of two parts. In the first (time-independent) part 
it is shown that the wave function of the discrete pseudo Laplacian has plane wave 
asymptotics, provided that the c-functions determining the orthogonality measure 
of the polynomials satisy certain analyticity requirements (guaranteeing that the 
spectrum of the discrete pseudo Laplacian is absolutely continuous) . This part of 
the discussion hinges on previous results describing the large-degree asymptotics of 
the class of multivariate orthogonal polynomials under consideration J33j I)4c . The 
second (time-dependent) part consists of a stationary phase analysis that permits 
proving the existence and unitarity of the wave operators and scattering operators 
given the plane wave asymptotics of the wave functions. Key ingredient of this part 
of the discussion is a stationary phase estimate from |RS[ p. 38-39] that controls the 
decay for t ±cxd of certain oscillatory integrals describing the difference between 
interacting and freely evolving wave packets. 

The paper is organized as follows. Section|21describes the construction of orthog- 
onal polynomials related to root systems. In Section 13 we introduce a commuting 
system of discrete pseudo Laplacians on the Weyl chamber diagonalized by the 
orthogonal polynomials in question. The wave operators and scattering operators 
for our discrete pseudo Laplacians are determined in Section ^ The stationary 
phase analysis that lies at the basis of the computation of these wave- and scat- 
tering operators is relegated to Sectional Finally, in Section El we specialize to the 
case of Macdonald polynomials and detail the scattering theory of the associated 
hyperbolic Ruijsenaars-Schneider type lattice Calogero-Moser models. Some key 
properties of the Macdonald polynomials invoked in Section have been collected 
in Appendix El at the end of the paper. For the reader's convenience, we have also 
included an index of notations in Appendix IbI 

Let us conclude this introduction by providing a brief description of what the 
main results amount to in the elementary (classical) situation of a root system of 
rank 1. Let c{z) be a zero- free analytic function on the disc \z\ < g, with g > \^ 
that is real- valued for z real and normalized such that c(0) = 1. We associate to 
c{z) an orthonormal basis of trigonometric polynomials Po{£,) t PiiO t P'^iC) i ■ ■ ■ for 
the Hilbert space L^((0,7r), that is obtained by applying the Gram- 

Schmidt process to the Fourier-cosine basis 1, cos(^), cos(2^), .... It is an immediate 
consequence of the three-term recurrence relation for the orthonormal polynomials 
Pf(^) that the wave function 



(1.1) 
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satisfies an eigenvalue equation of tlic form L'i' = 2cos(^)'I', where L represents a 
discrete (self-adjoint) Laplacian acting on lattice functions <j) E £'^{N) as 

L(t)e ^ ae(t)i+i + be(l)i + ae_i(t)i-i (0_i = 0), (1.2) 

with ai, be denoting the coefficients of the three-term recurrence relation. For c{z) = 
1, the polynomials Pe{^) amount to the Chebyshev polynomials of the second kind 
Ue{cos£,) = sin(£ -I- l)^/sin^, whence the wave function in Eq. (|l.l|l reduces in 
this case to the Fourier-sine kernel '^f\0 = 2sin(£ + 1)^. The Laplacian L ((T^ 
then amounts to a free Laplacian L*^"^ whose action on lattice functions is given by 

Theorem 14.11 (below') now states that for ^ ^ oo the wave function ^'f(^) (|1.1|) 
converges exponentially fast in L^((0, tt), (27r)^^d^) to the anti-symmetric combi- 
nation of plane waves 

*-(e) ^ si/2(^)e»(£+i)« _ s-i/2(^)e-»(^+i)«, (1.3) 

with s(0 = c(e-*«)/c(e*«). 

Furthermore, let us denote by F : ^ £^((0, tt), (27r)-idO and : 

f{K) ^ L^{{0,n),{2TT)-^d£,) the Fourier pairings with kernel and 4'f^(0, 

respectively: 

< ^^^TT , < ^^^TT . (1.4) 

Then Theorem 14. 21 and Corollary (below) state that the wave operators 0± = 
s — limt^ioo e**-^e~'*^*''' and the scattering operator S = fi^^il_ exist in and 
are given explicitly by the unitary operators ^± = o S^^l"^ o and S = 
(JF(°))"^ o5o where S denotes a unitary scattering matrix that is characterized 
by its multiplicative action on a wave packet (j) e L2((0,7r), {2Tr)-^d^) of the form 
{S4>){0 = si^OHO for < ^ < TT (with s(^) as defined just below Eq. ((T3|) l. 

Acknowledgments. Thanks are due to S.N.M. Ruijsenaars for several helpful 
discussions and to the referees for suggesting some improvements in the presenta- 
tion. 

2. Orthogonal Polynomials Related to Root Systems 

In this section we introduce a class of multivariate orthogonal polynomials related 
to root systems. For basic facts on root systems we refer to the standard works 

2.1. Polynomials on the Weyl Alcove. Let E, (•,•) be a real A^-dimensional 
Euclidean vector space and let il C E denote an irreducible crystallographic root 
system spanning E. We write Q and Q+ for the root lattice and its nonnegative 
semigroup generated by the positive roots 



Q = Spanz (R), Q+ = Spanp, (fi+ ) , 



(2.1) 
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and we write V and for the weight lattice its nonnegative cone of dominant 
weights 

V = {A e E I (A,a'') G Z, Va G i?}, (2.2a) 

V+ = {A G E I (A,a^) G N, Va G fi+}, (2.2b) 

where we have introduced the coroot = 2a/ {a, a). The algebra of (trigonomet- 
ric) polynomials on the Wcyl alcove 

A {e G E I < (e, a) < 27r, Va G JS+} (2.3) 

is spanned by the basis of the monomial symmetric functions 

where W C GL(E) denotes the Weyl group of the root system i?, = w(^), and 

\W\\ stands for the order of the stabilizer subgroup Wx = {w G | w(A) = A}. 

2.2. Factorized Weight Functions. We will now introduce a class of smooth 
weight functions on the Weyl alcove A that factorize over the root system R. To 
this end we write R^) = {a ^ R \ 2a ^ R} and i?i = {a G | f ^ R}. (So for 
a reduced root system one has that Rq = Ri = R and for the nonreduced root 
system R = BCn one has that Rq = Cjv and Ri — Sjv.) The weight functions 
under consideration are of the form 

A(^) = \ , (2.5a) 

with 

m= n ci.|(e-^<"'«>), (2.5b) 



where it assumed that the c- functions C|o,|(z) building C{C) (|2.5bp depend only on 
the length of the root a (so C|q,| (z) = c\p\ {z) if a and (3 lie on the same Weyl-orbit). 
For technical reasons, we will furthermore assume that these c-functions c\a\{z) 
are (i) analytic and zero- free on a closed disc Dg = {z G C | |z| < p} of radius 
Q > 1, (ii) normalized such that C|q|(0) — 1, and (iii ) real-valued for z G M (so 

2.3. Gram-Schmidt Orthogonalization. The technical conditions on the c-func- 
tions ensure that A(f) (|2.5a(l . Ij2.5b|l defines a smooth positive weight function on 
A (which extends analytically to a Weyl-group invariant function on E). We em- 
ploy this weight function to endow the space of trigonometric polynomials on the 
Weyl alcove with an inner product structure via embedding in the Hilbert space 
i2(A,A|^|2dC): 

(/,5)a = 4(A) J^fiOmMO \mfd^, yf,9 e L'{A,A\SfdO, 

(2.6) 

where g{£^) stands for the complex conjugate of .g(^), Vol(A) = /a > ^(0 
denotes the Weyl denominator 

S{0= n (e^<"-«>/2_e-^<".«)/2). (2.7) 
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Let ^ be a (partial) order of the dominant weights refining the dominance 
partial order 

fi ^ \- ^(zQ+ (2.8) 

such that the highest-weight spaces Span{r7ip}^£-p+^^_;;^, remain finite-dimensional 
for all A S V~^. By applying the Gram-Schmidt process to the partially or- 
dered monomial basis {mx}\^-p+ , we construct a normalized basis {Px}xe'P+ of 
L^(A, A|(5p di^) given by trigonometric polynomials of the form 

PxiO^ XeV+, (2.9a) 

with coefficients a^/i G C such that 

(where a\x > by convention). The Gram-Schmidt process guarantees that the 
polynomials Px, A S P"*" are orthogonal when comparable in the (partial) order 
>: (i.e. {Px,Pfj.)^ = when A >- ^ or A ^). Hence, a sufficient condition 
to ensure that our polynomials form an orthonormal basis of the Hilbert space 
L^(A, A|(5p da;) is to require the refinement >r of the dominance order ^ to be a 
linear ordering of . (The fact that the polynomials in Eqs. (|2.9a(l . I)2.9b|l form a 
complete set in L^(A, AjJp dx) is a consequence of the Stone- Weierstrass theorem.) 
In general, the orthonormal basis in question depends on the choice of such linear 
refinement. It will turn out below, however, that for our principal applications the 
c-functions are such that the orthogonality is already guaranteed when taking for 
>: simply the dominance ordering ^ (|2.8|) itself (in other words, in such case the 
construction results to be independent of the choice of the linear refinement). From 
now on we will always assume that we have fixed a sufficiently fine (partial) ordering 
>: so as to guarantee that the basis {Px}xe'P+ be orthogonal (i.e. (Px,Pf_i)^ ~ 
when X ^ i-i). 

2.4. Weyl Characters. The simplest example of the above construction is the 
special case with unit c-functions, i.e., with C|c|(z) = 1, Va G . The weight 
function then becomes of the form A(^) = 1 and the Gram-Schmidt process turns 
out to be independent of the choice of the refinement >i of ^ (i.e. in this case we 
may take >^ to be equal to ^ without restriction). The corresponding orthonormal 



polynomials PxiS.) amount to the celebrated Weyl characters jM8l IM4j 

P,(0=XA(e)^ri(e) ^(-l)-e^<''+^-^-\ XeV+, (2.10) 

where (—1)™ = det(u') and p = ^ J2aGR+ ^■ 

For later use, it will actually be convenient to extend the definition of the Weyl 
characters Xa(0 ™ Eq. (|2.1U|) to the case of nondominant weights A. It is immediate 
from this definition that for A G P \ V'^ 

x.(o = "i^-^!^;' (2.11) 

[0 if \Wp+x\ > 1, 

where, for /i G "P regular, Wfj_ € W denotes the unique Weyl group element such 
that w^ip.) e V+. 
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3. Discrete (Pseudo) Laplacians on the Weyl Chamber 

In this section we associate a commuting system of discrete pseudo Laplacians 
on to our orthonormal polynomials P\{^). 

3.1. Fourier Transform. Let H be the Hilbert space P{V^) of square-summable 
functions over the dominant cone equipped with the standard inner product 

if,9(^l'{'P+)), (3.1a) 

and let H be the Hilbert space L^(A,d^) of square-integrable functions over the 
Weyl alcove equipped with the normalized inner product 

if^9)n - |VF|Vol(A) A^'^^^^^^ (/,5ei'(A,dO). (3.1b) 
By construction, the functions 

^^{0 = ^'/'(mOPxiO, Aep+ (3.2) 

form an orthonormal basis of H. As a result, the mapping !F : H H given by 
4>x ^ 4>iO with 

= (0,*(O)« (3.3a) 

xev+ 

constitutes a unitary Hilbert space isomorphism between Ti. and Ti. The inverse 
mapping : H i-^ H takes the form (f){£^) — > (/)x with 

<j>x = {$,'^x)^ (3.3b) 

^ ^ <^(e)*A(e)de 



\W\Yol{A) J J, 

We will refer to J- as the Fourier transform associated to the polynomials Pa (C) • Iii 
the simplest case with unit c-functions, the wave functions amount to plane waves 
(cf. Section ITIjl 

*([')(^) ^ ^ (-l)™e^(p+^'«™). (3.4) 

The corresponding Fourier transform reduces to the conventional Fourier transform 
■.n^n of the form (t>x ^ <j>{0 with 

kO = E (-1)" E '^Ae-'^<^+^'«™>. (3.5a) 
The inverse transform (jF*^"')^^ :'H^ Ti'is then given by (/)(^) — > with 
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3.2. Pseudo Laplacians. To the basis of fundamental weights uji, . . . ,ujn gener- 
ating , we associate bounded multipUcation operators Ei, . . . ,i?jv in 'H of the 
form 

MO^ E exp(«('^'0), r^l,...,N, (3.6) 

where the sum is over ah weights in the Weyl orbit of cUr- The puUbacks of 
El, . . . ,En with respect to the Fourier transform define an integrable system 
of bounded commuting operators in H 

Lr = o E,.oT, r = l,...,N. (3.7) 

We will refer to the commutative algebra M[-Li, . . . , Ljv] generated by these oper- 
ators as the (algebra of) discrete pseudo Laplacians associated to the polynomials 
P\(^). It is immediate from its construction as the puUback of a multiplication 
operator in H (cf. Eq. H3.7|l l that the pseudo Laplacian Lr has a purely absolutely 
continuous spectrum in H given by the compact set cr{Lr) — {Er{^) \ ^ € A} C C. 
By acting with both sides of the operator equality L^T'^ — T^-^Er on (the complex 
conjugate of) an arbitrary element G 7i, we get 

i,(*A,0)7i = (^r^'A,^),^, V^eTi. (3.8a) 

In other words, the functions 'I'a(0 form a complete (as J- : Ti ^ H is a Hilbert 
space isomorphism) set of generalized joint eigenf unctions of our pseudo Laplacians, 
i.e. formally 

Lr^xiO = ^.(0*a(0- (3.8b) 
Here ^ G A plays the role of the spectral parameter and the weight A e is 
interpreted as the discrete geometric variable (i.e. the position variable). 

Note. i. Below we will sometimes write formal equalities of the form in Eq. (|3.8b|) 
that admit a rigorous interpretation of the form in Eq. I|3.8a|l upon taking the inner 
product (smearing) with an arbitrary (stationary) wave packet (/) G Ti.. 

Note. ii. In general the Laplacian Lr is not self-adjoint. Indeed, the adjoint L* is 
given by Ls with lOs ~ ~'WQ{LUr), where wq denotes the longest element of the Weyl 
group W (i.e., the unique Weyl group element wq such that wa{A) ~ —A). Thus 
Lr is self-adjoint if and only if wolcur) = —tOr- 

3.3. Localization. Let (p '■ 'P^ ^ C be square-summable a lattice function. The 
action of Lr on is of the form 

Lr(f'x = E ax^-r4>ti., (3.9) 

for certain coefficients a\^-^r G C. We will now show that in fact only a finite number 
of these coefficients is nonzero. 

Proposition 3.1 (Localization). The action of the pseudo Laplacian Lr on (f) E Ti. 

is of the form 

Lr4'X = X! °->^Kr4'tii CL\tJ.;r G C, (3.10a) 

where 

Plr = {a* G 7^+ I A* ^ A + cj^ and ^ - wo{ujr) h A}. (3.10b) 
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Proof. From the triangularity of the monomial expansion of Px{S,) it is immediate 
that 

The orthonormality furthermore implies that 
with LOs — —WQ^Wr) (cf. Note ii. above). Hence 

Since Er{0 = m^AO and *a(C) = A1/2(^)^(^)Pa(C), we conclude that 

Taking the innerproduct with an arbitrary wave packet (j) d Ti and comparison with 
the eigenvalue equation in Eq. H3.8a|l entails that 

i.e. formally (cf. Note i. above) 

The proposition then follows by the completeness of the generalized eigenfunctions 
^a(C); C a in the Hilbert space Ti (i.e. by the fact that the Fourier transform 
J- (|3.3afl . Ij3.3b|) constitutes a unitary Hilbert space isomorphism between Ti and 

n). □ 

A priori the cardinality of the set may be unbounded as a function of 

A G V'^ . Hence, in general our pseudo Laplacians need not be difference operators. 
If the ordering of the dominant weights > coincides with the dominance order ^, 
however, then it follows from the definition in Eq. (|2.8|l that the size of the set V^.^ 
is bounded by the number of weights in the interval {u ^ V \ W(){LOr) ^ v ^ w^}. 
Consequently, in this situation our pseudo Laplacian Lr is actually a difference 
operator in 7i. (We will refer in such case to Lr as a discrete Laplacian as opposed 
to merely a pseudo Laplacian.) 

Proposition 3.2 (Discrete Laplacians). When our ordering >: coincides with the 
dominance ordering ^ (I2.8|l . then the pseudo Laplacians in . . . , Lj^] are dis- 

crete difference operators in 7i. 

In the case of unit c- functions (cf. Section lT^ . our discrete Laplacians ii, . . . , Ljv 
amount to conventional free Laplacians L^^\ . . . , l''^^ over the dominant cone 7^+. 

Proposition 3.3 (Free Laplacians). //c|q|(z) = 1, Va G , then our discrete 
Laplacians Lr reduce to the free Laplacians 
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with the boundary condition that for fj, V \ 

(-l)'^^+-0»,+,(p+p)-P ^f \Wp+^\ = 1, 
z/ \Wp+^\ > 1 

( where Wp+^ denotes the Weyl permutation taking the regular weight p + ii to the 
dominant cone). 

Proof. As pointed out in Section 12.41 the case of unit c-functions corresponds to 
orthonormal polynomials -Pa(C) given by the Weyl characters X\{0- It is immediate 
from the explicit expression for x\ in El- l|2-10|l that the Weyl characters satisfy 
the well-known recurrence relations 

Starting from these recurrence relations, the proposition readily follows by repeating 
the arguments in the proof of Proposition 13. II The boundary condition stems from 
the property H2.11|l of the Weyl characters. □ 

By Proposition the eigenvalue equations in Eq. (|3.8b|l take the form 

aA;.;r*p(C)=4(0*A(C), l,...,^, (3.11a) 



E 



or equivalently 



E 



ax^^MO ^ EriOPxiO, l,...,iV (3.11b) 



(upon dividing out the trivial overall normalization factor 5{£,)/ \J C(^)C(— ^) on 
both sides). The latter equations admit an alternative interpretation as a system 
of recurrence relations (or Fieri formulas) for the polynomials P\{C)- 

Note. The above construction of the discrete (pseudo) Laplacians has its origin 
in the works of Macdonald [Mil IMS IMll- Specifically, for c\a\{z) = (1 - t^^^z) 
with —1 < < 1 the polynomials P\{£,) (|2.9a|) . Ij2.9bp amount to (the parameter 
deformations of) Macdonald's zonal spherical functions on p-adic Lie groups [Mil 
IM3| . The algebra of discrete Laplacians ]R[Li, . . . ,Lm\ corresponds in this case to 
the iC-spherical Hecke algebra of the p-adic Lie group. When C|q|(z) is given by a 
g-shifted factorial (cf. Eq. below), then the polynomials P\{C) specialize to 

the Macdonald polynomials (M2l IM3I IM4| . The discrete Laplacians appear in this 
context in Cherednik's double affine Hecke algebra as "coordinate multiplication 
operators" dual to Macdonald's difference operators |C2I IM4| . 



4. Time-Dependent Scattering Theory 

In this section we determine the wave operators and scattering operator asso- 
ciated to our discrete pseudo Laplacians. For background literature on scattering 
theory the reader is referred to e.g. Refs. |RS|Ip1It] . 
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4.1. Plane Wave Asymptotics. The dominant Weyl chamber is given by the 
open convex cone 

C+ = {x e E I (x,a) > 0, Va G (4.1) 

We will now describe the asymptotics of the wave function v1/a(C) H3.2|I diagonalizing 
the pseudo Laplacians Li, . . . , Ln (|3.7|I for A deep in the Weyl chamber, i.e., for A 
growing to infinity in such a way that (A, a^) +oo for all positive roots a G . 
To this end we define for A G 

m(A) EE min (A,a^). (4.2) 

aeR+ 

In previous work, it was shown that the strong L^-asymptotics of the polynomials 
PxlO for to(A) ^ oo is given by |R4llD3llD4j 

P^iO - E (-l)"'C(e^)e*<^+^'«™>. (4.3) 

More precisely, one has that 

II^A - Pr IIa - 0(e-"-(^)) as to(A) ^, (4.4) 

1/2 

where \\-\\^ = ("i")^ ^^'^ e > denotes a decay rate that depends on the radius g > 
1 of the analyticity disc of the c-functions C|q|(z) (see the technical assumptions 
in Section 

The idea of the proof in jR4l ID3I ID4j of this exponential convergence goes along 
the following lines. Firstly, a direct (constant term) computation reveals that 

(4.5a) 

I i IL fj, = A. 

Next, we denote by p(™(^)'(^) the polynomial approximation of the asymptotic 
function P^{£,) obtained by replacing the overall c- function C(^) by its Taylor 
polynomial of degree m(A). Then a combinatorial analysis shows that this poly- 
nomial approximation expands triangularly on the basis of monomial symmetric 
functions 

Pf(^»(0=.m,(a+ E ^Vm^a (4.5b) 

(for certain coefficients b\f^ G C). Moreover, the analyticity requirements on the 
c-functions C|a|(z) guarantee that 

P^{^) = P7(^)(e) + 0(e-™(^)) (4.5c) 

(this is because the technical conditions ensure that the Taylor coefficients of the 
c-function C|q|(z) decay exponentially fast). From Eqs. 14.5a|l - 14.5c|l one concludes 
that — up to an 0(e^""*^'^'') error term — the asymptotic function P^{S,) amounts 
to a monic polynomial obtained by performing the Gram-Schmidt process on the 
monomial symmetric basis with respect to the inner product (•, Oa- other words, 
the asymptotic functions coincide up to exponentially decaying error terms with the 
monic versions of the polynomials P\{£,) defined in Eqs. H2.9a|l . (|2.9b|l . The conver- 
gence in Eq. I|4.4|l now follows from the fact that the orthonormalized polynomials 
P\{£,) are asymptotically monic: a\x = 1 + 0(e~^™('^'). (This estimate for the 
leading coefficient in the monomial expansion of -Pa(0 follows starting from the 
equality a\\ — {P\, Px°)^, upon substituting H4.5c|l and expanding the polynomial 
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part P^^^\0 ill terms of the orthonormalized polynomials Pf^iO, /i ^ A, taking 
into account the orthogonality Ij2.9b|) .'l 

The asymptotic estimate in Eq. (|4.4|) for the polynomials -Pa(C) immediately 
gives rise to the following plane wave asymptotics for the wave functions ^'a(C) 

^TiO = A^/'imOPnO (4.6a) 



where 



with 



SUO = (4.6c) 

n ^vi(("'0) n wmt, (4.6d) 

^h((«,0)= 1 ' , (4.6e) 



(so S|a|((a,0) = S|a|(^("'0) = ^|a|(("'0) ^ud |s|a| ((«, 0) I = 1 ). 

Theorem 4.1 (Plane Wave Asymptotics). The wave function 5" a tends to the 
plane waves for X deep in the Weyl chamber: 

||*A-*rik = 0(e-""(^)) a,sm(A)^oo, 

w/iere || . ||^ = (t)^''^. 

We see from Theorem 14 . II that the asymptotics of the wave functions ^'a(0 for 
A deep in the Weyl chamber is given by an anti-symmetric combination of plane 
waves e*^'^'^^ with phase-factors that factorize over the root system in terms of 
one-dimensional c-functions. 

4.2. Scattering and Wave Operators. For any real multiplication operator 
E{£.) C R[Ei{C),...,En{0], let L = T-^oEoT&nd let L^o) = (J^W)-! oiJo jc-(o). 
In other words, the operators L C . . . , Lat] and L^'^^ C . . . , L^^\ are 

self-adjoint (pseudo) Laplacians in 71 such that (formally) 

L^xiO^m-^m and L(«)vl/(")(0 = ^(0*i°^(0. (4.7) 
(So the spectrum of L and L^^^ in H is absolutely continuous and given by the 
compact interval cr(L) = a{L^°^) = {E{£,) \ ^ G A}.) We wiU now describe the 
scattering of the interacting dynamics generated by the discrete pseudo Laplacian 
L with respect to the free dynamics generated by the discrete Laplacian L^^^ Let 
us to this end define the regular sector of the Weyl alcove as 

A„g = {e e A I {VE, a) / 0, Va e R+}. (4.8) 

Due to the analyticity of E{£,), the regular sector Arcg is an open dense subset of 
the Weyl alcove A. For every ^ G Arcg, there exists now a unique Weyl group 
element € W such that w^C^E) lies in the dominant Weyl chamber C"*" 1)4. 
Clearly, the Weyl-group valued function ^ ^ is constant on the connected 
components of Areg (by continuity). We are now in the position to define the 
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unitary multiplication operator Sl '■ 'H ^ Ti (the so-called scattering matrix) via 
its restriction to the dense subspace of (say) smooth complex test functions with 
compact support in Aj-cg: 

{SLm) = S^.iOm (<^eCo°°(A,eg)), (4.9) 

where Sw{C) is given byEq. j46d| . 

The main result of this paper is the following explicit formula for the wave 
operators and the scattering operator in terms of the scattering matrix Sl (|4.9|) 
and the Fourier transforms T (|3.3a(l . (|3.3bp and T^^^ (|3.5a(l . (|3.5bp . thus relating 
the long-time asymptotics of interacting dynamics e**^ to that of the free dynamics 
e'*^'"'. The proof, which is relegated to Section below, consists of a stationary 
phase analysis based on the asymptotic formula for the wave functions in Theorem 

Theorem 4.2 (Wave Operators). The operator limits 

n± = s- lim e'*^e-'*^*°' 

1 /2 

converge in the strong \\ • W^u-norm topology (where \\ • Wn — A ^'^^ 

corresponding wave operators : Ti Ti are given by the unitary operators 

Corollary 4.3 (Scattering Operator). The scattering operator Sl '■ H i-^ Ti, for 

the self-adjoint discrete pseudo Laplacian L G . . . , Ljsi] is given by the unitary 

operator 

5L^r!;^r!_ = (^("))-io5LO^("). 

We see from Corollarv l4.3l and Eqs. 14.6d|) . 1)4. 9|l that the scattering matrix for 
the self-adjoint discrete pseudo Laplacian L factorizes over the root system R. For 
the type A root systems. Theorem 14.21 and Corollarv 14.31 reproduce the results of 
Ruijsenaars in Ref. jR4j. 



5. Stationary Phase Analysis 

In this section the fundamental formulas for the wave operators stated in Theo- 
rem ^21 ^-re proven. To this end we employ a stationary phase method that gener- 
alizes Ruijsenaars' approach in Ref. R4 from the type A root systems to the case 
of arbitrary crystallographic root systems. Throughout this section the notational 
conventions of Sections |31 and 01 are adopted. 

5.1. Asymptotics of Wave Packets. Let us introduce the free wave packet 4>'^^^ [t) 
and the interacting wave packets (f)± (t) of the form 

0(0) (i) = (.F(°))-ie-^*^0, (5.1a) 
0±(i) = e-'*^ S]^^^^ (5.1b) 
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or more explicitly 

0fW = 4(A) 1^^-^^"' /^e'(^+^^^->--^(^)0(OdC, (5.2a) 

^.At) - ^^^^/^*^(Oe-*^«)>Sf/^(0^(e)de, (5.2b) 

with (j) € C^(Ai-cg). The following lemma states that the long-time asymptotics 
of the interacting wave packets (t>+(t) and (f'-it) for t +oo and t — > — oo, re- 
spectively, coincides with the corresponding asymptotics of the free wave packet 
0(0) (t). 

Proposition 5.1 (Asymptotic Freedom). Fort ±oo, the difference between the 
interacting wave packet (j)± (t) and the free wave packet (^("^ (t) tends to zero: 

Vk>0: \\(t)±{t)-<j>^'^\t)\\n^O{l/\t\'') ast^±oo. 

Before proving Proposition 15.11 (cf. below) , let us first infer that Theorem 14.21 
arises as an immediate consequence. Indeed, since the space of test functions 
C^{Arcg) is dense in H and the operators in question are unitary, to validate The- 
orem 14.21 it is sufficient to demonstrate that for = with cj) S C^(Ai.eg) 

lim ||e'*^e-^*^''V-f^±'^llw = 0, 

t — >±QO 

where Q± = o Sf^^"^ o JTC^). From the unitarity of e**'^ and the intertwining 
relations 

-itL^"^ (''Zr(0)\-1 (--27(0) \-l -its „ J -itL t--1 t--1 -itE 

e o [j-^ ') — l^j-^ ') o e and e o y- — j- o e , 

it is clear that 

= U^'Ht) - ^±it)\u 

whence the theorem follows from Proposition l5.ll 

5.2. The Classical Wave Packet. To prove Proposition 15.11 we may assume 
without restricting generality that the test function (p has in fact compact sup- 
port inside a connected component of A^cg. In this situation there thus exists a 
unique Weyl group element w ^ W such that w(Vij) lies inside the dominant 
Weyl chamber C+ (|4.1|) for all ^ in the support of cj). Let Vdas C E be an open 
bounded neighborhood of the compact range of classical wave-packet velocities 
Ran^(V£^) = {VE{£,) \ £, € Supp(0)} staying away from the walls of the Weyl 
chamber w~^{C'^) in the sense that there exists a lower-bound £ > such that 
(CiiijO;^) > e for all C G Vdas and all a S . We will now introduce a classical 
wave packet that is finitely supported on the following i-dependent region of the 
cone of dominant weights V'^ 

p+ f{Ae7'+ |p+Aetl«(Velas)} for i > 0, 

ciasl } ]^{XeP+\p + Xe twow(V,i^)} for i < 0. ^ ' ' 

Because of dimensional considerations, it is clear that the cardinality of the support 
^cias(^) grows at most polynomially in t 

lnla.(0l=O(i^) forl^Hoo. (5.4) 
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The classical wave packet is defined as 



(5.5) 



|l^|Vol(A) 



${0d^ for A G V+Jt) and t > 0, 



|M/|Vol(A) Jj^ 



e«{p+A,c„„^)-»t£;(0^(^)dC for A e V+i^^{t) and t < 0, 
otherwise. 



The next lemma compares the long-time asymptotics of the free wave packet 
c/)^^'' (t) with that of the classical wave packet 

Lemma 5.2. For t — > ±00, the difference between the free wave packet (j)^^^{t) and 
the classical wave packet tends to zero: 

Vk>0: ||0(°)(t)-0("'"'')(i)||w = O(l/|i|'^) ast-^ ±00. 

Proof. It is immediate from the definitions of the wave packets under consideration 
that 



^(0)(,)„^(clas)(^) 



1 



\W\Yo\{A) 



EM) 



s'^^+^'^-^-'^^^^^COdC, (5.6a) 



with 



W 



W\{iD} ifAeP+^,(t)andi>0, 
W \ {wqw} if A G V+^^{t) and t < 0, 
W otherwise. 



(5.6b) 



The proof of the lemma now hinges on a stationary phase estimate extracted from 
the Corollary of Theorem XI. 14 in Ref. ESI p. 38-39], which states that for any 
fc > there exists a (positive) constant Ck such that 



„j{x,e)-jt£;(c) 



< 



Ck 



(l + |x| + |<| 



for all X e E and t £ R such that 



X ^ iVclas- 



(5.7a) 



(5.7b) 



Indeed, invoking of the stationary phase estimate in Eqs. (|5.7a|) . Ij5.7b|l with k > 
N/2 and x = w~^{p + A), reveals that the norm of the difference between the 
wave packets given by Eqs. (|5.6a|l . H5.6b|l in the Hilbert space H is 0(l/|t|'^~^/^) 
as i ^ ±00. (Notice in this connection that w^^{p + A) G iVdas if and only if 



A G 7'+^,(i) and weW\W.) 



□ 



5.3. The Asymptotic Wave Packet. Let us define asymptotic wave packets 
4>^^\t) that are obtained from the interacting wave packets 4i±{t) upon replac- 
ing the Fourier kernel 4'a(C) by its plane wave asymptotics 5'^°°''(C) 



1 



iy|Voi(A) 



*i°°^(Oe-**^«^5f^'(O0(Od^, (5.8a) 
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or more explicitly 

The next lemnia states that the long-time behavior of the asymptotic wave packets 
is governed by the classical wave packet (t) H5.5(l . 

Lemma 5.3. For t —>■ ±00, the difference between the asymptotic wave packet 
(t)j- (t) and the classical wave packet tends to zero: 

Vk>0: = 0(l/|ir) as t^ ±00. 

Proof. The proof of Lemma l5.2l aDDlies verbatim, upon replacing (t) by (j)^^\t) 
and the introduction of minor modifications in the formulas so as to incorporate 
the additional (smooth) factors C{^w)/C{iw) and C{$w)/C{^wow), respectively. □ 

Let p^'^^^'^'' : 7-^ denote the orthogonal projection onto the finite-dimensional 
subspace fiKi.sit)) ^ 

It is clear from the definition of the classical wave packet that p/'^^^'^ (^('^i^^) (i)) ^ 
0(cias) j-^-j^ As a consequence, we get from from Lemma 15.31 upon projection onto 
^'(^ls(0) that 

Vk>0: ||Pt^''"'V±°^(t) = 0(l/|ir) asi-.±oo. (5.10) 

Our final lemma states that the long-time asymptotics of the interacting wave 
packet 4>±{t) coincides with that of the asymptotic wave packet <l)'^\t). 

Lemma 5.4. For t ±00, the difference between the interacting wave packet 
4>±{i) and the asymptotic wave packet (j)^^\t) tends to zero: 

Vk>0: U±{t)-<t)^^\t)\\H^O{\/\t\^) as t^ ±00. 
Proof. From the definitions it is immediate that 

0±,A(t) - 0S(t) = {e'^'^Sf/'l^, *f ^)^. 

Hence 

||pf-)(</,±(t)-4°°)(t))||?, = \{e-''^Sr/^,^,-¥r')n\' 

(by the Cauchy-Schwarz inequality). Now, since |'Pcias(^)l — 0{t'^) and m(A) > 
\t\e~ 1 for A e '^cias(Oj '^6 conclude from this estimate combined with Theorem l4.1l 
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that llp^''^^'*^-' (i) — {t)j \\-H converges to zero exponentially fast for t — > ±00, 
so in particular 

V«;>0: llP/"'"'' (</)±(t)- 4°°' (t)) |lw = 0(1/10 ast^ioo. (5.11) 

The lemma now follows from the vanishing of the tails (Id — pl-'^^^^'')(j)[^\t) and 
(Id - Pt^''''^)(f>±{t) for t -> ±00: 

||(Id-Pf'"'Vi°°'(t)||w = 0(l/|tr) ast^icx), (5.12a) 

||(Id-P/'''''^)0±(t)||w 0(l/|t|'') asi->±oo. (5.12b) 

Notice in this connection that the limiting relation in Eq. (|5.12a(l is immediate from 
Lemma lOl and Eq. (|5.10(l . and that the limiting relation in Eq. I)5.12b|l follows by 
compairing the norm equality 

with the norm estimate for t — > ±00 

\\pt'''''U±mn ^''^ iipf^^VL^^(oii«+o(i/itr) 
E..m ||^(-^)||^ + o(i/|tr) 

LemmaO || (t) || ^ + 0( 1/ |i | - ) 

which entails that 



||(Id-P/^'^^V±Wlh = \^U±itWn - \\4'''''U±(tWn = 0{l/\tr^'). 



□ 



5.4. Proof of Proposition l5.ll After these preparations, the proof of Proposition 
15. II reduces to the telescope 

||0±(t)-0(")(i)||« (5.13) 
and the application of Lemmas 15 . 21 15 . 31 and 15 .41 



6. Application: Scattering of Hyperbolic Lattice Calogero-Moser 

Models 

In this section we specialize our c-functions so as to describe the scattering of the 
hyperbolic Ruijsenaars-Schneider type lattice Calogero-Moser models associated 
with the Macdonald polynomials. Initially, viz. in the first three subsections, it 
will be assumed that our root system R be reduced (so Rq = Ri = R) except 
when explicitly stated otherwise. In the fourth subsection, we then indicate how 
the results extend to the case of a nonreduced root system (so R = BCn, Rq = Cn 
and ill = Bn)- We end our study of the lattice Calogero-Moser models in the fifth 
subsection by providing some illuminating additional details describing what the 
results boil down to in the simplest situation of a root system of rank one. 
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6.1. Macdonald Wave Function. For c- functions of tlie form 

c\a\[z) = — , g = e , 5|„|,s>0, (6.1) 

witii {z;q)oo = Y[^=oi^ ~ tlie weight function A(^) (|2.5a|l ~ (j2.5b|) becomes 

(Tire positivity restrictions on the parameters g\a\ and s guarantee that the c- 
function H6.1|l meets the technical requirements stated in Section |21) Our polyno- 
mials P\{Oj ^ £ 'P^ ^^'^^ amount to the orthonormalized Macdonald polynomials 
[M2limilM4] 



^a(C) = ^Ai/2(A)P,(0, (6.3a) 



where 



A^^^ C+{Pg)C jPg) _ C jPg) 

- TTTl 1 —TV' -^0 - 77-7 T, (6.c!b) 

C+{Pg + X)C [Pg+X) C+{pg) 



with Pg ^ ^ J2aeR+ 9\a\ a and 



C±(x) = n c±|((x,a^)), (6.3c) 

Here Pa(0 denotes the Macdonald polynomial 

Pa (6 = caPa(0, ca - ^^^f;^' (6-4a) 

where 

Pa(0="^a(0+ E ca^™m(0, (6.4b) 

with coefhcients ca/^ G C such that 

(pa, "^/i)A =0 for ^ ^ A. (6.4c) 

For the Macdonald weight function A(^) (|6.2|) . the coefficients c\fj, turn out to 
vanish when A and p are not comparable in the dominance ordering |M3| . In other 
words, in this case one may take the ordering >^ to be equal to the dominance 
order ^ (|2.8(l without restricting generality. Explicit formulas for the expansion 
coefficients ca^ when V ^ Q (i.e. excluding the root systems £^8, and G2) can 
be found in Ref. [PLMj . 

We thus arrive at the following formula for the wave function ^a(C) H3.2|) in 
terms of Macdonald polynomials. 

Proposition 6.1 (Macdonald Wave Function). For R reduced and c-functions 
given by C|q|(z) 1)6.1(1 . the wave function ^'a(0 I|3.2(I reads explicitly 



a(0 = -^AV2(A)AV2(e)^(e)P,(C), 
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R 


minuscule 


quasi 


A AT ■ 


w\ , • • • , UN 


U\ 


Bn ■ 


UN 


^i, 


Cn '■ 


Ui 




Dn : 


Wi, a;jv_i, UN 


W2, 


Ee : 


Ul.Ud 


W2, 


Er : 






Es : 




Ws, 


F4 : 




W4, 


G2 : 




^1, 


BCn ■ 







i-minuscule 



Table 1. Minuscule and Quasi- Minuscule Weights. 



with Pa(C) denoting the Macdonald polynomial characterized by Eqs. (|6.4a|| - (|6.4c|l . 

6.2. Macdonald-Ruijsenaars Laplacian. Let us recall that a nonzero weight 
TT G is called minuscule if (tt, a^} < 1 for all a G R^ and that it is called 
quasi-minus cule if (tt, a^) < 1 for all a G /J^ \ {tt} (and it is not minuscule). The 
number of minuscule weights is equal to the index of Q in 'P minus 1 (so for Eg, F4 
and G2 there are none). As regards to the quasi- minuscule weights: there is always 
just one and it is given by tt = ao, where ag is the maximal root of the dual root 
system R^ = {a^ | a e R}. For the readers convenience, we have included a list 
of the (quasi-)minuscule weights for each root system in Table ^ (where we have 
adopted the standard numbering of the fundamental weights in accordance with 
Refs. lEUHni). 

To a (quasi-)minuscule weight tt we associate the multiplication operator i?^ : 
H H given by 

EAO= E exp(z(i.,0). (6.5) 

iyeW{-!r)UW{--iT) 

The Macdonald-Ruijsenaars Laplacian is now defined as the puUback : 7i > 7i 
of i?7r with respect to the Fourier transform J- 

= oE^oT. (6.6) 

By ProDOsition lH.2l the Macdonald-Ruijsenaars Laplacian Ltt G ffi[ii, . . • , Ln] con- 
stitutes a difference operator in Ti. The following proposition provides its explicit 
action on lattice functions over the dominant cone . 

Proposition 6.2 (Macdonald-Ruijsenaars Laplacian). For R reduced andu (quasi- 
)minuscule, the action of the Macdonald-Ruijsenaars Laplacian L^^ on a ( square- 
summable) lattice function (j) : C is given by 

L^(p\ = E^{p^g) + 

i'GW(7r)UW'(-7r) 
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whe 



Y ( \ _ TT (g|a| + (x,a'^) : sinhs)(^.av) 
"^""^ " ((x,av) :sinh,)<,,„v) 

CKt-Tt 

(!^,a'')>0 

„V\ 



n 



sinhf (g|a| + (x, a^)) 
sinhf((x,aV)) ^ 



n sinh|(g|a| + (x, a^)) sinh§(l + (x, g^)) 

sinhf ((x,aV)) sinh f (1 + (x, a^^) ' 

{v,a')=2 

£^7r(x) = ^ exp(s(i/,x}), 

i/eVK(7r)UlV(-7r) 

wit/i (z : sinhs)m = Jlfco^ sinh §(z + £) and Pg^^ J2aeR+ 9\a\<^ 



V 



Proof. It is a straightforward consequence of the definitions that the c-functions 
satisiy the difference equations 

Ch(2; + 1)_ sinh(f) cp^|(a; + l) _ sinhf(l + x-5i„|) 



cj^|(a:) sinhf (g|„| +a;)' sinhf(l + a;) 

With the aid of these difference equations it is not difficult to verify the fundamental 
functional relation 

A(x + J^)Vl^(pg+x + j/) = A(x)K(Pg+x), i^eW{T:). (6.7) 

From the recurrence relation for the Macdonald polynomials exhibited in Eq. ljA.6a|) 
of the Appendix, it is now readily inferred — upon invoking the functional relation 
(|6.7(l specialized to x A with A and A + dominant — that the Macdonald wave 
function 'I'aIO lj^3SJ-lj^3iI) satisfies the eigenvalue equation 

E + >^)vli'{pg + A + ,y)^x+AO - VAP9 + A)*a(0) 

uew{7T) 

i^Giy(7r) 

Combining with the corresponding eigenvalue equation in which tt is replaced by 
—wo{'k) (where wq is the longest element of the Weyl group), leads us to the eigen- 
value equation for i^: 

E {V^'^Pa + m-J{P9 + A + ^)^x+AS.) - VuiPa + A)*a(0 



i/GVF(7r)UVK(-7r) 



(K(?)-K(pP)*a(0- 



The proposition now follows from the completeness of the Macdonald wave functions 
*a(C), C G A in the Hilbert space □ 
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When TT is minuscule we have that 



and that 



E 



K(x)-£;.(pp (6.8b) 

i'eW(7r)UW(-7r) 

(by the Macdonald identity in Eq. (|A.5p of the Appendix). As a consequence, the 
action of the Macdonald-Ruijsenaars Laplacian in Proposition l() . 21 simphfies in this 
situation to 

L^^x= J2 V,'^\pg + X)vli\pg + X + iy)q}x+,. (6.8c) 

ueW(T:)UW(-TT) 

For the root system Ajv, all fundamental weights are minuscule (cf. Table 
Hence, in this special case the discrete Laplacians Li, . . . , Lj^ of SectionOare given 
by the Macdonald-Ruijsenaars Laplacians Lt^ H6.8cf) with tt = uj, j = 1, . . . , A^. 
The operators in question correspond to the commuting quantum integrals of the 
hyperbolic relativistic lattice Calogero-Moser model due to Ruijsenaars pRll ilR3| . 
For the other root systems, only a small part of the discrete Laplacians Li, . . . , Ljv 
can be made explicit by means of the Macdonald-Ruijsenaars Laplacian of Proposi- 
tion l^^l^nd Tabled In principle, the higher-order commuting Laplacians may be 
constructed with the aid of the corresponding Dunkl-Chcrednik difference-reflection 
operators |C1IIC2| . however, at present explicit formulas for a set of generators for 
the algebra of commuting Laplacians . . . ,L]\[] are available only in the case 

of the classical root systems |D2I |5] . 

It is of course a consequence of our construction that the algebra of Laplacians 
M[-Li, . . . ,Ln] consists of bounded self-adjoint operators in the Hilbert space 7i. For 
the Macdonald-Ruijsenaars Laplacian i^, one can also check this fact independently 
directly from the explicit action in Propostion l6.2l 

Note. If TT is quasi-minuscule then — tt G W{tt) (as tt is a root). Thus, in this case 
W{—Tr) — W{tt). The same simplification also occurs for tt not necessarily quasi- 
minuscule when —1 e W (i.e. when the longest Weyl-group element Wq equals — 1). 
This is the case for the root systems B^f, Cjv, Dn (^ > 4, even), Ej, Eg, F^, G2 
and BCn, but it is not the case for the root systems Ajv (A^ > 2), Z)jv (-^ > 3, 
odd) and Eq. 

6.3. Scattering Matrix. When g^^i 1, Va G R, the Macdonald c-functions 
C|q| (z) (|6.1|I specialize to the unit c-function. The Macdonald-Ruijsenaars Laplacian 
of Proposition 16 . 21 reduces in this limit to the free Laplacian it'^'' — (jF^^^)"^ o 
E^ o jc-(o) given by 

u£W(-it)\JW(-tt) 

with boundary conditions as stipulated in Proposition 13 . 31 The following proposi- 
tion provides a somewhat more explicit characterization of these boundary condi- 
tions (in the case of tt (quasi-)minuscule). 
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Proposition 6.3 (Action of the Free Laplacian). For n (quasi-) minuscule, the 
action of the free Laplacian L^^^ is of the form 

with n„{X) — Q if t: is minuscule, and with nT;{\) denoting the number of short 
simple roots aj perpendicular to X if tt is quasi-minuscule (where, by convention, 
all roots qualify as short if R is simply laced). 

Proof. Starting point is tlic action of Li"' in Eq. H6.9() with boundary conditions 
as detailed in Proposition 13.31 li X + v ^ , then there exists a simple root 
aj such that (A + J^, aJ) < 0. Hence, since A e 7^+ and v e W{'k) U W{-it) 
with TT (quasi-)niinuscule, it follows that we are in either one of the following three 
situations: 

(i) (A,aJ) =0 and {v,a)) = -1, 

(ii) (A, a)) = and {v, a)) = -2, 

(iii) (A,aJ) = 1 and {v,a)) = -2. 

It is not difficult to see that in the first and last situation the weight p + A + i/ is 
stabilized by the simple reflection r^^ . Indeed, we get 

raj{p + X + u) = p + X ^- u - {p + X + V, aJ) = p + X + v 

(where we exploited the fact that (p, a^) = 1 for a simple). It thus follows that in 
these two cases the stabilizer of p + A + is nontrivial, whence the corresponding 
term (px+u in Eq- vanishes by the boundary condition in Proposition 13.31 

The second situation occurs only when tt is quasi-minuscule. Clearly we must 
then have that u = —aj, whence (p + A + i',a'j) = —1 and (p + A + J^, a^) = 
1 + (A, aj^) — {aj,a^) > for k ^ j (where we exploited the fact that {a, (3^) < 
for a, /3 simple and distinct). It thus follows that the weight p + A + is regular 
and that the Weyl permutation Wp+\+i, taking it to the dominant cone is given by 
the simple reflection r^^. . Indeed, we now get 

raj{p + X + i/)=p + X + u-{p + X + iy, aj) aj — p + X. 

Invoking of the boundary condition in Proposition 13 . 31 then reveals that the corre- 
sponding term 0A+iy in Eq. (|6.9(l is equal to Now, every simple root aj in the 
Weyl orbit of tt for which (A, aJ) = gives rise to such a contribution (px-aj = —px 
in the action on the r.h.s. ofEq. H().9|) . Furthermore, since a quasi-minuscule weight 
TT is a short root of R (as tt — ag with q:q denoting the maximal root of R , whence 
is long and ag is short), it is clear that the nonzero contributions in question 
occur precisely at all simple short roots perpendicular to A. □ 

Our main application of the scattering formalism in Sectionals the following ex- 
plicit formula for the scattering and wave operators for the lattice Calogero-Moser 
system, relating the long-time asymptotics of the dynamics of the Macdonald- 
Ruijsenaars Laplacian to that of the free Laplacian -LI-^'. 

Theorem 6.4 (Lattice Calogero-Moser Scattering for R Reduced). The wave op- 
erators il± s — limj^ioo e'*^''e~'*^'' ' and the scattering operator Sl^ = ft^^fl- 
for the Macdonald-Ruijsenaars Laplacian L^^ in relation to the free Laplacian L^^ 
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are of the form stated in Theorem \4--°A O'f^d Corollary with a unitary scattering 
matrix Sl^{£,) given by Eqs. (14. 9|) . (|4.6d|) and 

For the type A root systems Theorem l6.4l is due to Ruijsenaars jR4) . The scat- 
tering of the corresponding classical-mechanical system was analyzed previously in 
Ref. 1^21. 



Note. The asymptotics of the Macdonald wave function ^'a(C) in Proposition lti.ll is 
governed by Theorem 14. II fand Eqs. (|4.t)a|) - H4.6e|l ') with a scattering matrix taken 
from Theorem l6.4l 

6.4. Extension to Nonreduced Root Systems. We will now indicate how the 
results of Subsections I6.1H6.3I should be adapted so as to include the case of a 
nonreduced root system (viz. R = BCn, Ro — Cn, Ri = Bj^ and W amounts 
to the hyperoctahedral group Sm k Z^). In short, the bottom line is that all 
results carry over to the case of nonreduced root systems upon passing from the 
Macdonald polynomials to the Macdonald-Koornwinder multivariate Askey- Wilson 
polynomials [KllU2l l5]. More specifically, by picking c-functions C|q|(z), a G Ri of 
the form 

{qSz;q)oo „ , 

— ^ — ^ — tor a long, 

— — ^ — for a short, 

[qz^;q)oa 

where g = e"-* and s, 5, go, ■ • ■ ,ff3 > (and with (zi, . . . , Zk;q)oo = (2i;g)oo • • • {zk;q)oo), 
we end up with a weight function A(^) H2.5al) - I|2.5b|l given by 

{qae}i"^i);q) 



^(^) = n 7^^?l% (6-11) 



aeRi 
a long 



a short 

The polynomials P\{^), A G amount in this case to orthonormalized Macdonald- 
Koornwinder polynomials K, D2, S . The polynomials in question are again of the 
form in Eqs. (|6.3a|l . H6.3b|l and Eqs. I|6.4a|l - H6.4c|l . but now with 

C±(x)= n c±|((x,a)), (6.12a) 

c+|(x)= (6.12b) 
f^W2(9!l!i9k for a long, 

(<? ;9)oo 

qi9o+gi+g2+g3)x/2 ^ 

(n9o+x _„gi+x g2 + l/2+x g3+l/2+x \ 

for a short, 

(r^;9)oo 
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(6.12c) 



ox/2 (g '^;g)oo 

q{ga+gi+g2+g3)x/2 ^ 



for a long, 



(ql-So+a:^ —gl-Sl+a:^ gl/2-ff2+a;^ _ql/2-g3+x.^ 



for a short, 



where we have distinguished dual parameters 5, • ■ • i53 that are related to the 
parameters g, ^q, . . . ,^3 via the linear relations 



(6.13a) 



1 ] 






1 


1 


1 






( 9.\ 


91 


1 


1 


1 


-1 


-1 






9i 


92 


~ 2 


1 


-1 


1 


-1 






92 


V 93 ) 




I 1 


-1 


-1 


1 


) 




\ 93 J 



9 = 9. 

and with the vectors pg and now taken to be 

P9 = f H " + 50 X! 



a long 



a<£RJ 
a short 



Pg 



^ Oi + go ^ a. 



(6.13b) 
(6.13c) 



a long 



a short 



We thus arrive at the following formula for the wave function 5'a(0 J3.2|I in terms 
of Macdonald-Koornwindcr polynomials. 

Proposition 6.5 (Macdonald-Koornwinder Wave Function). For R nonreduced 
and c- functions given by c\a\(z) l|6.1()|l . the wave function 5'a(0 I|3.2|1 reads explicitly 



= ^Ai/2(A)Ai/2(e)5(^)p,(e), 



with Pa(C) denoting the Macdonald-Koornwinder polynomial characterized by Eqs. 
(|?najl -(|Oc | . 

From the second-order recurrence relation for the Macdonald-Koornwinder poly- 
nomials |D2| , we now obtain the following formula for the action of the Macdonald- 
Koornwinder Laplacian = J-~^ o o T , associated to the first fundamental 
weight TT = cji (which is a quasi- minuscule weight for R — BCn, cf. Table 

Proposition 6.6 (Macdonald-Koornwinder Laplacian). For R nonreduced and 
TT = oji, the action of the Macdonald-Koor winder Laplacian on a (square- 
summable) lattice function (j) : —>■ C is given by 

E {^'^^(Pa + ^)V-i\pg + A + iy)^x+. - V,{pg + A)0a), 
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wh- 



ere 



K(x) = n !^^Mk_M^x 



a long^ (i^,a)— 1 



sinh §((x, a)) 



-r-r sinh I (ffo + (x, q) ) cosh § (gi + (x, a) ) 
smh§((x,a}) coshf((x,a}) 

a short, (i/,q) — 1 

sinh §(ff2 + ^ + (x, g}) cosh §(33 + ^ + (x, a)) 
sinh |(i + (x, a)) cosh |(i + (x, a)) 

EA^) = X! cxp(s(j/,x)), 

and TOt/i pg and given by Eqs. (|6.13b|) and (|6.13c|l . respectively. 

For 5 — > 1 andjo, ... ,33 — > 1/2 (so g,go — > 1 and 31,32,33 — > 0), the c- 
functions C|q|(z) (|6.1U|I tend to 1 (recaU in this connection the duphcation formula 

{z'^;q)ac, = {z, —z,q^^^z, —q-^/^z;q)ao)- The Macdonald-Koornwinder Laplacian 
then reduces to the free Laplacian 

4°)0A= E ^^+- (6.14) 

Application of the scattering formalism of Section^ now produces the following 
scattering and wave operators relating the long-time asymptotics of the dynamics 
of the Macdonald-Koornwinder Laplacian to that of the free Laplacian L?-* 
(EH. 

Theorem 6.7 (Lattice Calogero-Moser Scattering for R Nonreduced). The wave 
operators ^± = s— lim(^±oo e**^''e~'*^'' ' and the scattering operator Sl^ = fJ^^fi- 
for the Macdonald-Koornwinder Laplacian in relation to the free Laplacian L^^^ 
(|6.14|l are of the form stated in Theorem \4-'^ and Corollary \4-^ with a unitary 
scattering matrix Sl^ given by Eqs. H4.9|) . 1)4. 6dl) and 

■S|a|((",0) = 

r (9e^<"'«>;g)oo (g^e-^<«'«>;g)oo 



(g9e'<"^«);5)oo (ge-'<"'«);g)oo 

(ge2'<"'«);q). 



for a long, 



for a short. 



Note. The asymptotics of the Macdonald-Koornwinder wave function 5'a(C) in 
Proposition 16.51 is governed by Theorem 14.11 fand Eqs. (|4.t)a|) - (|4.6e|l ') with scat- 
tering matrices taken from Theorem 16. 71 

6.5. Example: The Rank-One Case. It is quite instructive to exhibit the results 
of this section in somewhat further detail for simplest case of a root system of rank 
one. We will restrict attention the case of a nonreduced root system (i.e. BCi), 
since the reduced case (viz. Ai) can be recovered from it via a specialization 
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of the parameters (corresponding to a standard reduction from the Askey- Wilson 
polynomials to the g-ultraspherical polynomials |AWIIGS] ). 

In this situation the Macdonald-Koornwinder wave function takes the explicit 
basic hypergeometric form 

^iiO = ^^'^'ie)A'^'imO'PeiO, ^eN, ee(0,7r), (6.15) 



where 



^fo = (6.16a) 

m = "^^'t'[''^ (6.16b) 

Mo = WTTT' (6-16c) 

SiO = 2sin(0, (6.16d) 



with 



^^^^ = (^^^^^ ' ('-^'^^ 

c+{x) = qi9o+9i+g2+a,)x/2 ^ (6.17b) 

(^q90+x ^ —q9i+^^ q92 + l/'i+x ^ _qg3 + l/2+x. q-j^ 

('Z^^;9)oo 

C-{x) = q(90+9l+92+g3)x/2 ^ (6.17c) 

(g^+^";g)oo 

(^ql-ga+x^ _ql-gi+x^ ql/2-g2+x^ _ql/2-g3+x. q^^ ' 

and with Pe{£,) denoting the Askey- Wilson polynomial |AWI IGb^ 

/ q^^ ^q^9o+^ ^qUo^ii ^qgog-ii, \ 

P£(0 = 4$3 y _qgo+g\ ^ ^3o+S2 + l/2^ 1^90+93 + 1/2 5 9, 9 j ■ (6.18) 

Here we have employed standard notation from the theory of basic hypergeometric 
series (tS] 

/ ai,...,as \ (fli, . . . ,as;g)„ 



with (a;g)„ = Ilfc^oll " ^-^d (ai, . . .,as;q)n = {ai;q)n ■ ■ ■ {as\q)n- 

The asymptotics of the wave function 5';(C) for I — > cxo is given by (cf. also 

vI/-(^) = si/2(0e'('+i)« - ri/2(^)e-'(^+i)€, (6.19a) 

with 

^(0 = (6.19b) 

The free plane waves 

*£°'(0 EH boil in this case down to the Fourier sine 

kernel 

(0 = 2 sin(£ + 1% ^ G N, e e (0, tt). (6.20) 
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The corresponding Fourier pairings T : ^^(N) ^ L^((0,7r), (27r) ^d^) and T'^^^ : 
i^((0,7r), (27r)^^d^) together with their inversion formulas are given by 



1 

2^ 



^sN- (6.21a) 







and 




(6.21b) 



10 



respectively (where we have omitted the complex conjugations because the relevant 
kernel functions '^li^) and £^re real- valued as a consequence of the fact that 

-1 e w^Zs). 

The Macdonald-Koornwinder Laplacian L = T^^ o Eo T and the free Laplacian 
lio) ^ (jr(o))-i o£;o jc-(o) associated to the multiplication operator E{^) = 2cos(f) 
act on on lattice functions <j) £ respectively as 

Lq^e = V'/^go+£)V'/^i~go'~i-l)(l)i+i+ (6.22a) 
V'^^i-go - e)V'/^{go + I - \)<i>i-i + 
(2cosh(s5o) - ■(^(ffo +^) - (1 - k,o)y{-9o ~ t))<t'i: 



with 



= sinhf(go + .)cosh|(g,+.) 
smh|(a;) cosh | (a;) 

sinh f (52 + 5 + a;) cosh f (53 + 5 + 2:) 
sinh|(i+a;) cosh|(i+a;) 



and as 



= 0^+1 +0^-1, (6.23) 
with the boundary condition =0. 

The specialization of Theorem 16.71 to the case = 1 now states that the wave 
operators Sl± = s — hmt_±oo e^'^e"**-^*"' and the scattering operator Sl = il^^fi- 
exist in ^■^(N), and are moreover of the form il± = JF^^ o 5^^^^ o and Sl = 
(^_^(o)^-i oSl o respectively, with Sl being a unitary scattering matrix whose 
multiplicative action on a wave packet (f) S i^((0, tt), (27r)~^d^) is given by 

{SLm)^'^m for0<e<7r. (6.24) 
c(?) 

Appendix A. Properties of the Macdonald Polynomials 

This appendix serves to list a number of key properties of the Macdonald poly- 
nomials Pa(0) ^ defined by Eqs. H6.4a|l - (|6.4c|) . We used these properties in 
Section 13 to build the Macdonald wave function and to the determine the explicit 
action of the Macdonald-Ruijsenaars Laplacian. For proofs of the statements below 
and for further theory concerning the Macdonald polynomials the reader is referred 
to the seminal works of Macdonald and Cherednik |M2l |M3l |M4| [UTl (see also 
[C] for a different approach). Throughout this appendix it is assumed that our 
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root system R be reduced. For the extension of the statements below to the case of 
nonreduced loot systems the reader is referred to Refs. |M3I IM4l IkI ID2I lOkl IS] . 

The Macdonald polynomials Pa(C) l|t)-4a|) - (|6.4c|l are normalized such that they 
satisfy the Specialization Formula 

Pa(*spP ^ 1. (A.l) 
In this normalization the Orthogonality Relations read 

The specialization formula in Eq. ljA.l|) amounts to the special case /i = of the 
more general Symmetry Relation 

P^i^^iP^, + A*)) = Pf i^^iPa + A)), (A.3) 

where Pf (C) and P^'' (^) refer to the Macdonald polynomials associated to the root 
system R and the the dual root system i?^ , respectively (so A and ^ are dominant 
weights of R and R^ , respectively). 

For any (quasi-) minuscule weight tt of i?^, we have a corresponding Macdonald 
Difference Equation given by 

T-r («s.g|a| + : sin)(^,„) /e , ■ n r. f \ a \ 

^ n ((e,a):sin),.., (P.(g + z..) - P.(0) (A.4a) 

{y,a)>0 

= E ('Z^'^''^"^^ Pa (0, 

where (z : sin)™ = Ilfco^ ^{z + is£). If the weight tt is minuscule (so [(tt, a)! < 1, 
Va S R), then this difference equation simplifies to 

i^eW{7r) aeR 2V\';'"/; ueW(,7r) 

(z^,q) — 1 

because of the Macdonald Identity (for tt minuscule) 

{u,a) — l 

Combination of the Macdonald difference equation in Eq. ljA.4all and the sym- 
metry relation in Eq. ljA.3|) leads to the Recurrence Relation (or Fieri Formula) 

(e'(--«> _g(-^pr))P^(^) = (A.6a) 
.ewM aeR ((Ps + A,av) :smh,)<.,„v) 
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where tt is now a (quasi-)minuscule weight of R (and (z : sinhs), 
£)))■ In the minuscule case this recurrence relation reduces to 




sinh(§(z+ 



sinh§(ff|a| + {pg + A,a^)) 
sinh|((pg + A,Q!V)) 



Pa+.(0- (A.6b) 



Appendix B. Index of Notations 



This Appendix provides a list of notations ordered according to the sections in 
which they were first introduced. 



SectionO E, (•, •) , R, R+ , Q, Q+, V, V+ , a', A, 77.^(0, W, W^, \Wx\. 

Sectional Ho, -Ri, A(e), ^(0, C|„|(z). 

SectionO (•,.)a, Vo1(A), (5(0, ^, Pa(0, «Am- 

SectionO XA(e), (-1)"', P, ^^m- 

SectionO (•,•)«, ^, (•,-)7i, *A(e), ^, (0, 

SectionO w^, Er{Cj, W{-), Lr, cr{Lr), wq- 

Section axfj,;r, V^.^, 

SectionO C+, m(A), P^{i), \\-\\^, Pf^\0, ^TiO, S^iO, 5|a|((a,0), IMI^^ 

SectionO EiO, L^, a{L), A,eg, w^, Sl, \\ ■ \\n, Sl- 

SectionO (t)^°\t), (j)±{t). 

SectionO Vcias, V+^,{t), 4,(-'^\t), W. 

SectionO -P/''"^- 

SectionO 5|a|, {z;q)oo, A(A), A/'o, C±(x), pg, c^|(x), PA(e). 

SectionO ^, a^, -R^, ^.(0, i^, K(x), E^{^), {z : sinh,)™, p^. 

Section 1^31 L^^\ n-ir(A), aj, r^j, Sl^, Sl^- 

Section|^2 g, gr, (zi, . . . , z^; g) 

oo, 3: 9r- 

SectionESl s+i^s, {a;q)n, (ai, . . . , a^; g)„. 



AS. M.J. Ablowitz and H. Segur, Solitons and the Inverse Scattering Transform, SIAM Stud- 
ies in Applied Mathematics 4, Society for Industrial and Applied Mathematics (SIAM), 
Philadelphia, Pa., 1981. 

AW. R. Askey and J. Wilson, Some basic hypergeometric orthogonal polynomials that generalize 

Jacobi polynomials, Mem. Amer. Math. Soc. 54 (1985), No. 319. 
B. N. Bourbaki, Groupes et algebres de Lie, Chapitres 4^6, Hermann, Paris, 1968. 



References 



30 



J.F. VAN DIEJEN 



C. O.A. Chalykh, Macdonald polynomials and algebraic intcgrability, Adv. Math. 166, 2002, 
193-259. 

CI. I. Clicrednik, Double afBne Hccke algebras and Macdonald's conjectures, Ann. Math. 141 
(1995), 191-216. 

C2. , Macdonald's evaluation conjectures and difference Fourier transform. Invent. Math. 

122 (1995), 119-145. 

Dl. J.F. van Diejen, Integrability of difference Calogero-Moser systems, J. Math. Phys. 35 
(1994), 2983-3004. 

D2. , Self-dual Koornwinder-Macdonald polynomials, Invent. Math. 126 (1996), 319- 

339. 

D3. , Asymptotic analysis of (partially) orthogonal polynomials associated with root 

systems, Internat. Math. Res. Notices 2003 (7), 387-410. 
D4. , Asymptotics of multivariate orthogonal polynomials with hyperoctahedral symme- 
try, in: Jack, Hall- Littlewood and Macdonald Polynomials (V.B. Kuznetsov and S. Sahi, 

eds.), Contemp. Math., Amer. Math. Soc, Providence, RI, (to appear). 
DLM. J.F. van Diejen, L. Lapointe, and J. Morse, Determinantal construction of orthogonal 

polynomials associated with root systems, Compositio Math. 140 (2004), 255-273. 
DV. J.F. van Diejen and L. Vinet (eds.), Calogero-Moser-Sutherland Models, CRM Series in 

Mathematical Physics, Springer- Verlag, New York, 2000. 
FT. L.D. Faddcev and L.A. Takhtajan, Harniltonian Methods in the Theory of Solitons, Springer 

Series in Soviet Mathematics, Springer- Verlag, Berlin, 1987. 
OS. G. Gasper and M. Rahman, Basic Hypergeometric Series, Encyclopedia of Mathematics 

and Its Applications 34, Cambridge University Press, Cambridge, 1990. 
Ga. M. Gaudin, La Fonction d'Onde de Bethe, Masson, Paris, 1983. 
Gu. E. Gutkin, Integrable systems with delta-potential, Duke Math. J. 49 (1982), 1-21. 
HO. G.J. Heckman and E.M. Opdam, Yang's system of particles and Hecke algebras, Ann. Math. 

145 (1997), 139-173; erratum ibid. 146 (1997), 749-750. 
HS. G. Heckman and H. Schlichtkrull, Harmonic Analysis and Special Functions on Symmetric 

Spaces, Perspectives in Mathematics 16, Academic Press, Inc., San Diego, CA, 1994. 
Hu. J.E. Humphreys, Introduction to Lie Algebras and Representation Theory, Springer- Verlag, 

New York, 1972. 

I. M.E.H. Ismail, Asymptotics of the Askey- Wilson and y-Jacobi jjolynomials, SIAM J. Math. 

Anal. 17 (1986), 1475-1482. 
IW. M.E.H. Ismail and J. A. Wilson, Asymptotic and generating relations for the g-Jacobi and 

4*3 polynomials, J. Approx. Theory 36 (1982), 43-54. 
K. T.H. Koornwinder, Askey- Wilson polynomials for root systems of type EC, in: Hypergeo- 
metric Functions on Domains of Positivity, Jack Polynomials, and Applications (D. St. P. 

Richards, ed.), Contemp. Math. 138, Amer. Math. Soc., Providence, RI, 1992, pp. 189-204. 
KBI. V.E. Korepin, N.M. Bogoliubov, and A.G. Izergin, Quantum Inverse Scattering Method 

and Correlation Functions, Cambridge University Press, Cambridge, 1993. 
Ku. P.P. Kulish, Factorization of the classical and quantum S matrix and conservation laws, 

Theor. Math. Phys. 26 (1976), 132-137. 
Ml. I.G. Macdonald, Spherical Functions of p-Adic Type, Publ. of the Ramanujan Inst., No. 2, 

1971. 



M2. , Symmetric Functions and Hall Polynomials, Clarendon Press, Oxford, 1995. 

M3. , Orthogonal polynomials associated with root systems, Sem. Lothar. Combin. 45 

(2000/01), Art. B45a, 40 pp. (electronic). 

M4. , Affine Hecke Algebras and Orthogonal Polynomials, Cambridge U.P., Cambridge, 

2003. 

M. D.C. Mattis (ed.). The Many-Body Problem: An Encyclopedia of Exactly Solved Models in 

One Dimension, World Scientific, Singapore, 1994. 
N. A.G. Newell, Solitons in Mathematics and Physics, CBMS-NSF Regional Conference Se- 



ries in Applied Mathematics 48, Society for Industrial and Applied Mathematics (SIAM), 
Philadelphia, PA, 1985. 

N-Z. S. Novikov, S.V. Manakov, L.P. Pitaevskii, and V.E. Zakharov, Theory of Solitons. The In- 
verse Scattering Method, Contemporary Soviet Mathematics, Consultants Bureau [Plenum] , 
New York, 1984. 



SCATTERING ON A WEYL CHAMBER 



31 



Ok. A. Okounkov, BC-typc interpolation Macdonald polynomials and binomial formula for 
Koornwindcr polynomials. Transform. Groups 3 (1998), 181-207. 

OP. M. A. Olslianctsky and A. M. Perelomov, Quantum intcgrablc systems related to Lie alge- 
bras, Phys. Rep. 94 (1983), 313-404. 

O. E.M. Opdam, Lecture Notes on Dunkl Operators for Real and Complex Reflection Groups, 
MSJ Memoirs 8, Mathematical Society of Japan, Tokyo, 2000. 

Ox. S. Oxford, The Hamiltonian of the Quantized Nonlinear Schrddinger Equation, Ph. D. 
Thesis, UCLA, 1979. 

P. D.B. Pearson, Quantum Scattering and Spectral Theory, Academic Press, London, 1988. 
RS. M. Reed and B. Simon, Methods of Modern Mathematical Physics. III. Scattering Theory, 

Academic Press, New York-London, 1979. 
RSc. S.N.M. Ruijsenaars and H. Schneider, A new class of intcgrablc systems and its relation to 

solitons, Ann. Phys. (NY) 170 (1986), 370-405. 
Rl. S.N.M. Ruijsenaars, Finite-dimensional soliton systems, in: Integrable and Superintegrable 

Systems (B. Kupershmidt, ed.). World Sci. Publishing, Teaneck, NJ, 1990, pp. 165-206. 
R2. , Action-angle maps and scattering theory for some finite-dimensional integrable 

systems. IIL Sutherland type systems and their duals, Publ. Res. Inst. Math. Sci. 31 (1995), 

247-353. 

R3. , Systems of Calogero-Moser type, in: Particles and Fields (G. Semenoff and L. 

Vinet, eds.), CRM Series in Mathematical Physics, Springer- Verlag, New York, 1999, pp. 

251-352. 

R4. , Factorized weight functions vs. factorized scattering. Comm. Math. Phys. 228 

(2002), 467-494. 

S. S. Sahi, Nonsymmetric Koornwinder polynomials and duality, Ann. Math. 150 (1999), 
267-282. 

SCM. A.C. Scott, F.Y.F. Chu, and D.W. McLaughlin, The soliton: a new concept in applied 

science, Proc. IEEE 61 (1973), 1443-1483. 
T. W. Thirring, A Course of Mathematical Physics, Vol. Ill: Quantum Mechanics of Atoms 

and Molecules, Springer- Verlag, New York, 1981. 



iNSTiTUTO DE Matematica y FfsiCA, Universidad de Talca, Casilla 747, Talca, Chile 



